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Bounding  superposed  on-off  sources  - 
Variability  ordering  and  majorization  to  the 

rescue 

Armand  M.  Makowski  * 


Abstract 

We  consider  the  problem  of  bounding  the  loss  rate  of  the  aggregation  of 
independent  on-off  sources  in  a  bufferless  model  by  the  loss  rate  resulting 
from  the  aggregation  of  i.i.d.  on-off  sources.  This  is  done  through  a  unified 
framework  based  on  the  interplay  of  well-known  results  from  the  theory  of 
variability  orderings  with  the  concept  of  majorization  ordering.  In  particular, 
we  use  a  basic  comparison  result  to  readily  derive  a  bound  of  Rasmussen  et 
al.  for  heterogeneous  sources  and  an  upper  bound  of  Mao  and  Habibi  for 
homogeneous  sources,  and  to  discuss  a  second  upper  bound  proposed  by 
these  authors.  It  is  argued  that  this  conjectured  upperbound  is  too  tight  in 
general,  and  should  be  replaced  by  new  and  provably  correct  upper  bound. 


1  Introduction 

Traffic  burstiness  has  long  been  considered  a  key  factor  for  provisioning  link  and 
buffer  resources  at  ATM  multiplexers.  In  a  first  step,  these  issues  can  be  addressed 
with  the  help  of  a  simple  bufferless  model  fed  by  fluid-like  input  traffic.  An  infor¬ 
mation  source  is  then  characterized  by  its  R+-valued  rate  process  { Bit),  t  >  0}, 
so  that  the  source  bursts  at  time  t  >  0  with  an  instantaneous  rate  of  Bit)  bps,  say 
for  sake  of  definiteness.  For  obvious  practical  reasons,  it  is  customary  to  require 
the  constraint  0  <  Bit)  <  P  (t  >  0)  where  P  is  the  peak  rate  of  the  source. 

*  Department  of  Electrical  and  Computer  Engineering,  and  the  Institute  for  Systems  Research, 
University  of  Maryland,  College  Park,  MD  20742.  E-mail:  armand@isr.umd.edu. 
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1.1  Loss  rates 

In  most  situations  of  interest,  the  rate  process  {/?(£),  t  >  0}  can  be  assumed 
ergodic  (as  we  do  from  now  on)  in  the  sense  that  for  all  x  >  0, 

lim  —  [  1  \R(t)  <  x]  dt  —  P  [R  <  x]  a.s.  (1) 

T— >oo  T  JO 

for  some  R+ -valued  rv  R.  If  the  rate  process  is  stationary  and  ergodic,  then  (1) 
holds  with  the  steady- state  rate  variable  R  determined  through  the  weak  conver¬ 
gence  R(t)  =ri  R.  Under  (1)  the  source  admits  an  average  rate  given  by 

m(R)  :  =  lim  —  /  R(t)dt  =  E  [R]  a.s.  (2) 

T->o o  T  Jo 

If  the  traffic  is  offered  for  transmission  over  a  a  link  operating  at  C  bps,  only 
min (R(t),C)  bps  can  be  accommodated,  and  in  the  absence  of  any  buffer,  the 
remaining  (R(t)  —  C)  '  1  bps  represents  the  instantaneous  loss  rate  over  that  link. 
Under  (1)  the  (average)  loss  rate  of  the  source  {R(t)  ,  t>  0}  over  the  C  bps  link 
is  well  defined  and  given  by 

L(fi;  C)  :=  lim  i  f  (R(i)  -  C)+dt 

r— >oo  l  Jo 

=  E[(i?-(7)+]  a.s.  (3) 

1.2  Multiplexing  sources 

While  the  definition  (3)  for  L (f?;  C)  might  appear  too  poor  a  marker  of  source 
behavior  to  be  of  any  use,  its  evaluation  is  nevertheless  helpful  either  for  di¬ 
mensioning  link  capacity  or  as  the  basis  for  a  Call  Admission  Control  (CAC) 
procedure  [1,  2,  7].  In  the  latter  instance,  traffic  carried  on  the  link  is  typically 
obtained  by  multiplexing  several  independent  information  sources.  If  N  sources 
{Rn(t),  t  >  0}  (n  =  1, . . . ,  N)  are  multiplexed  on  a  link  operating  at  C  bps,  the 
total  instantaneous  rate  is  then  given  by 

R(t)  =  Ri(t)  +  .  .  .  +  i?Af(f),  t  >  0. 

Under  appropriate  ergodic  assumptions,  it  follows  that 

L  (R-C)  =  L  {R1  +  ...  +  Rn-C) 

=  E  [(i?!  +  ...  +  Rn-  C)+]  (4) 

1  We  write  x+  =  max(a:,  0)  for  any  scalar  x. 
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where  the  mutually  independent  rvs  Ri, ... ,  Rjy  are  the  steady-state  rate  variables 
for  the  component  sources. 

As  indicated  already  in  [2,  3],  evaluating  (4)  can  be  computationally  pro¬ 
hibitive  even  in  the  simplest  of  cases  due  to  the  large  number  of  sources  that  need 
to  be  multiplexed  at  any  given  time.  This  difficulty  is  further  exacerbated  when 
the  component  sources  are  statistically  dissimilar  (as  is  the  case  in  practice)  [2]. 
This  state  of  affairs  has  prompted  a  search  for  upper  bounds  on  loss  rates  which 
are  computationally  efficient,  and  yet  sufficiently  tight  to  provide  good  approxi¬ 
mations. 

1.3  On-off  sources 

Most  of  these  efforts  have  been  carried  out  for  the  class  of  on-off  sources  (e.g.,  [1, 
2,  3,  7]).  A  source  with  rate  process  {R(t)  ,t>  0}  is  said  to  be  a  (generalized)  on- 
off  source  if  R(t)  alternates  between  two  states,  namely  R(t)  =  0  (resp.  R(t)  = 
P)  when  the  source  is  silent  (resp.  active)  at  time  t  >  0.  Under  the  ergodic 
assumption  (1),  such  an  on-off  source  admits  a  steady- state  rate  R  with  finite 
range  {0,  P}.  In  fact,  it  is  easy  to  see  that 

P[R  =  P]  =  1-  P[R  =  0]=  f(R) 

where  f(R)  is  the  activity  factor  of  the  source  defined  by 

f(R)  ■—  lim  —  f  1  \R(t)  >  0]  dt  a.s.  (5) 

T-*oo  T  JO 

For  on-off  sources,  we  have  m(R )  =  f(R)P,  so  that  such  sources  are  fully  (and 
equivalently)  characterized  by  either  of  the  pairs  (P,  f(R ))  or  ( P ,  m(R)).  We  find 
it  useful  to  represent  the  steady  state  rate  R  of  the  on-off  source  { Bit)  ,  t  >  0} 
with  peak  rate  P  and  activity  factor  f(R)  as 

R=st  PB(f(R ))2 

where  for  p  in  [0,1],  let  B(p)  denote  asn  (0, 1}- valued  (Bernoulli)  random  vari¬ 
able  (rv)  with  P  [B(p)  =  1]  —  p.  We  then  refer  to  such  an  on-off  source  as  the 
on-off  source  ( P ,  B(f(R))). 

1.4  Earlier  bounds  and  new  results 

Consider  A"  independent  on-off  sources  (Pn,  B(fn ))  with  peak  rate  Pn  and  activity 
factor  fn  (n  —  1, . . . ,  N)\  the  resulting  steady  state  rate  for  the  aggregate  traffic  is 

N 

Y.PnBn(Jn)  (6) 

_  n= 1 

2For  two  R-valued  rvs  X  and  Y  with  the  same  distribution,  we  write  X  =stY. 
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where  B\  (/, ), . . . ,  BN(fN)  are  independent  Bernoulli  rvs.  The  upper  bounds  de¬ 
rived  in  the  literature  on  loss  rates  (4)  for  the  aggregate  traffic  (6)  can  be  inter¬ 
preted  as  loss  rates  for  an  aggregation  of  fewer,  say  L  <  N,  i.i.d.  on-off  sources 
with  common  peak  rate  Pnew  and  activity  factor  /new.  The  resulting  steady  state 
rate  for  the  aggregate  traffic  is  now 

L 

EPnew^(/new)  (7) 

t=  1 

where  Bi(fnew), . . . ,  BL(fnew )  are  i.i.d.  Bernoulli  rvs.  As  will  become  apparent 
in  the  forthcoming  sections,  the  validity  of  the  comparison 

N  L 

L(E  PnBntfn)  :  C)  <  L(£  ^new  Wnew);  C),  C>  0  (8) 

n= 1  t=  1 

entails  tradeoffs  in  that  a  smaller  value  of  L  (desirable  for  obvious  computational 
reasons)  corresponds  to  a  larger  value  for  Pnew  (not  desirable  as  it  leads  to  looser 
bounds). 

Rasmussen  et  al.  [7,  p.  353]  conjectured  that  when  the  sources  in  (6)  have 
identical  peak  rates,  say  P,  but  possibly  different  activity  parameters,  the  aggre¬ 
gation  (7)  of  N  homogeneous  on-off  sources  with  identical  peak  rate  Pnew  =  P 
and  activity  parameter  /new  =  Ar~  1  (/,  +  . . .  +  fN ),  provides  an  upper  bound. 
This  conjecture  was  recently  established  by  Mao  and  Habibi  [3,  Thm.  1]  from 
basic  principles.  These  authors  also  establish  another  upper  bound  [3,  Thm.  3], 
this  time  for  N  homogeneous  on-off  sources,  by  replacing  them  with  a  reduced 
number  of  homogeneous  on-off  sources.  Finally,  they  conjecture  the  validity  of 
an  upper  bound  [3,  Conjecture  1]  which  generalizes  both  the  upper  bound  of  Ras¬ 
mussen  et  al.  and  their  upper  bound;  a  full  discussion  of  this  combined  upper 
bound  can  be  found  in  [4,  Thm.  3,  p.  127]. 

Here,  we  revisit  these  upper  bounds  by  establishing  a  general  comparison  re¬ 
sult  for  weighted  sums  of  independent  Bernoulli  rvs  [Proposition  4].  We  show 
how  this  general  result  provides  a  unified  vehicle  for  readily  deriving  the  bound 
of  Rasmussen  et  al.  for  heterogeneous  sources  and  the  upper  bound  of  Mao  and 
Habibi  for  homogeneous  sources,  and  for  discussing  their  second  upper  bound. 
It  is  argued  that  this  conjectured  upperbound  is  too  tight  in  general.  We  then  use 
Proposition  4  to  generate  a  new  upper  bound  which  is  provably  correct. 

The  proper  framework  for  addressing  these  issues  (and  similar  comparisons 
more  generally)  is  one  that  combines  stochastic  orderings  [9]  with  the  notion  of 
majorization  [6]:  The  variability  orderings  we  use  are  tailor-made  for  comparing 
loss  rates  as  in  (8),  while  majorization  is  useful  for  formally  comparing  degrees 
of  heterogeneity.  The  relevant  definitions  and  facts  are  given  in  Section  2.  This 
is  followed  in  Section  3  by  a  discussion  of  three  simple  operations  that  reduce 
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variability;  this  material  readily  yields  the  general  comparison  result  in  Section  4. 
Applications  of  the  general  result  are  presented  in  Section  5. 


2  Stochastic  orderings  and  majorization 

The  basic  tools  are  introduced  in  this  section. 

2.1  Variability  orderings 

For  IR-valued  rvs  X  and  Y,  we  say  that  X  is  smaller  than  Y  in  the  convex  (resp. 
increasing  convex)  ordering  if 

E  b(X)]  <  E  [p(y)l  (9) 

for  all  mappings  ip  :  R  — >  1R  which  are  convex  (resp.  increasing  and  convex) 
provided  the  expectations  in  (9)  exist;  we  write  X  <cx  Y  (resp.  X  <lcx  Y ).  We 
refer  to  these  orderings  as  the  variability  orderings.  Additional  material  on  these 
orderings  can  be  found  in  the  monographs  [8]  and  [9]. 

2.2  Key  facts 

We  now  present  well-known  facts  that  help  shape  the  approach  taken  here.  First, 
an  equivalent  definition  of  the  convex  increasing  ordering  [9,  Thm.  1.3.1,  p.  9]. 

Proposition  1  For  it -valued  rvs  X  and  Y  with  finite  expectations,  we  have  X  <lcx 
Y  if  and  only  if 


E 


(X 


<  E 


(Y 


a  e  1R. 


Proposition  1  makes  it  clear  why  the  variability  orderings  are  likely  vehicles  for 
carrying  out  the  comparisons  discussed  earlier.  Put  simply,  establishing  the  com¬ 
parison  L (Ri,C)  <  L (R2;  C )  for  all  values  of  C  between  the  loss  rates  of  two 
information  sources  with  steady-state  rates  Ri  and  R2  is  equivalent  to  the  com¬ 
parison  Ri  <icx  R2. 

Next,  we  explore  the  impact  of  the  constraint  E  [X]  =  E  \Y]  [9,  Thm.  1.3.1, 
p.  9], 

Proposition  2  ForH-valued  rvs  X  and  Y  with  finite  expectations,  wehaveX  <cx 
Y  if  and  only  if  X  <icx  Y  and  E  [X]  =  E  [Y]. 
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Finally,  the  convex  ordering  is  closed  under  independent  addition  [9,  p.  9]. 

Proposition  3  Consider  two  sets  of  mutually  independent  1R -valued  rvs  X1, . . .  ,XN 
and  Y\, ,  YN.  If  Xn  <cx  Yn  for  each  n  —  1, . . . ,  N,  then 

X\  +  . . .  +  Xn  Ycx  Y\  +  . . .  +  Yn- 


2.3  Majorization 


Let  K  denote  some  given  positive  integer.  For  any  vector  x  =  (.r  i , ....  xk) 
in  1Ra,  let  £(i)  <  X{2)  <  •  •  •  <  X[k)  denote  the  components  of  x  arranged  in 
increasing  order.  For  vectors  x  and  y  in  RA ,  we  say  that  x  is  majorized  by  y.  and 
write  x  -<  y,  whenever  the  conditions 

k  k 

J2x(i)  ^  k  =  l,2,...,K  (10) 

i— 1  i=l 


hold  with 


K  K 

Y,Xi  =  &' 


(ID 


i=l  i= 1 

Additional  information  regarding  majorization  can  be  found  in  the  monograph  [6] . 
Note  that  for  any  x  in  RA ,  we  have  xave  -<  x  with  e  =  (1, . . . ,  1)  in  RA ,  and 


xav  =  —  {xi  +  . . .  +  xK)- 

K 


3  Reducing  variability 

Below  we  identify  three  operations  that  reduce  variability,  thus  leading  to  com¬ 
parisons  in  the  ordering  <cx. 


3.1  Normalized  Bernoulli  rvs 

We  begin  with  a  comparison  result  for  renormalized  Bernoulli  rvs.  Recall  that  for 
p  in  [0, 1],  let  B(p)  denote  an  (0, 1}- valued  rv  with  P  [B(p)  =  1]  =  p. 

Lemma  1  The  collection  of  rvs  { p~lB{p ),  p  e  (0, 1]}  is  monotone  decreasing 
in  the  convex  ordering,  i.e., 

—  CxP  ~'B(p),  p  <  q. 
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In  other  words,  increasing  p  makes  p  1  B(jp)  less  variable. 
Proof.  We  need  to  show  that 


E 


<p(q  50?)) 


<  E 


<P(P  b(p )) 


p  <  q 


for  any  convex  mapping  ip  :  R  — >  R,  where 


(12) 


E 


‘ Pip  b(p )) 


p(<p(p  x)  -^(o))  +  ^(o),  pe(o,  l]. 


Hence,  it  suffices  to  establish  (12)  for  convex  mappings  <p  :  R  — >  R  such  that 
99(0)  =  0.  However,  under  this  constraint,  it  is  well  known  that  x  — >  x~l<p(x)  is 
non-decreasing  on  (0,  00)  and  the  conclusion  follows.  ■ 


3.2  Heterogeneity  decreases  variability 

For  p  in  [0, 1] A ,  we  define  the  rv  Sk(p )  as  the  sum 

Sk(p)  =  J2  BkiPk) 

k=  1 

where  the  Bernoulli  rvs  -Bi(pi),  •  •  ■ ,  bk{pk )  are  assumed  mutually  independent. 

Lemma  2  For  vectors  p  and  q  in  [0, 1]A,  it  holds  that  S^iq)  <cx  Sk{p )  when¬ 
ever  p  -<  q. 


Proof.  For  any  integer-convex  mapping  <p  :  IN  — >  R,  we  define  the  mapping 
:  [0, 1]A  — >  R  by 

*k(p)  =  E  [<p(SK(p))] ,  p  G  [0, if.  (13) 

It  is  well  known  [6,  F.l,  p.  360]  that  the  mapping  <I>  K  is  Schur-concave  in  that  the 
condition  p  -<  q  implies  ^k(q)  <  $k(p),  and  the  conclusion  Sx(q)  <cx  Sk ip) 
follows  from  the  definition  of  the  convex  ordering  <cx.  ■ 


The  next  result,  originally  due  to  Hoeffding  [6,  p.  359],  is  an  immediate  con¬ 
sequence  of  Lemma  2. 

Lemma  3  For  any  vector  p  in  [0, 1]A,  it  holds  that  Sk{p )  <cx  <SV(paVe)  where 
Pav  =  jfipi  +  •  •  -+Pk)- 
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3.3  Linear  combinations 

Let  {Xn,  n  —  1,2,.. .}  denote  a  sequence  of  i.i.d.  IR-valued  rvs.  The  following 
result  is  an  easy  consequence  of  Proposition  B.2  in  [6,  p.  287];  see  also  B.2.b  in 
[6,  p.  288]. 

Lemma  4  For  each  positive  integer  K,  it  holds  that 

K  I< 

'y  &kXk  5; cx  y  v  bfoXfo  (14) 

k= 1  fc=l 


whenever  a  -<  b  in  1R  . 

An  immediate  corollary  to  Lemma  4  is  obtained  by  taking  positive  integers 
L  <  K,  and  observing  that  a  -<  bin  [0,  1] h  with 

a  —  K~1(  1,  and  b  =  L_1(l, . . . ,  1,  0, . . . ,  0). 

Lemma  5  For  positive  integers  L  <  K,  it  holds  that 


1  k  1  l 

Xk  <cx  j  5Z  xt- 

n  k= i  L  i=  i 


(15) 


This  last  result  was  first  derived  by  Marshall  and  Proschan  [6,  B.2.c,  p.  288],  and 
formalizes  the  notion  that  averaging  decreases  variability. 


4  The  main  result 

Consider  N  independent  on-off  sources  as  described  in  Section  1.3,  where  for 
each  n  =  1, . . . ,  N,  the  nth  source  (Pn,  Bn(fn ))  has  peak  rate  Pn  and  activity 
factor  fn  so  that  its  average  rate  mn  is  given  by 

Win  Pnfn- 

As  these  N  sources  are  multiplexed,  the  resulting  total  average  rate  is  simply 

W, total  =  mi  +  ■■■  +  mN.  (16) 

Proposition  4  With  P*  selected  so  that 

max  Pn  :=  Pmax  <  P*,  (17) 

n=l,...,N 
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set 


r  := 


'^total 

iVP*  ' 


For  any  positive  integer  L  <  N,  it  holds  that 


N 


NP * 


EPnBn(fn)  <«r  E  W 


(18) 


(19) 


n=l  *=1 

where  the  rvs  Bi(f*), . . . ,  BL(f*)  are  i.i.d.  Bernoulli  rvs. 

Thus,  the  aggregation  of  heterogeneous  independent  on-off  sources  can  be 
upper  bounded  in  the  sense  of  the  convex  ordering  by  an  aggregation  of  fewer 
related  i.i.d  on-off  sources. 

Proof.  For  each  n  —  1, , . . ,  N,  define 

,*  Pnf 

■'  n  ‘  p *  -I  n 

and  note  that 

Pnfn  =  P*f ^  =  mn, 

so  that  /*  lies  in  (0, 1]  since  /*  <  fn.  From  this  last  equality  we  conclude  by 
Lemma  1  that 

<c,  /r‘Bn(/„*).  (20) 

where  we  take  the  Bernoulli  rvs  Bx  (/*), . . . ,  BN(f^)  to  be  mutually  independent 
rvs. 

With  this  in  mind,  we  now  get 


n= 1 


E  PnBn{fn)  =  E  Pnfn  Bn(  fn 

n= 1 
N 

y  py:  (k'bm 


n= 1 
N 


—  CX 


y  p*fn  (/r'B»(/, 

n=l 

N 

p *  £  B„(JD 


n= 1 


where  the  inequality  follows  from  (20)  via  Lemma  3. 
Next,  we  observe  that 


1  N  1  N  rn 

__  \  ’  f*  —  _  \  '  n  —  f* 

]\J  Jn  —  l\T  2^  p*  —  J 

iv  n= 1  iv  n=l 


(21) 
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Invoking  Lemma  3,  with  i.i.d.  Bernoulli  rvs  . . . ,  BN(f  *),  we  find  that 


N 

P'  E  BnU'n) 


n=  1 


<cx  P*£Pn(/*) 

n=  1 

=  NP’  (/  E ,«..(/')) 
<«  ^-p*(^Eb<(/*)) 


(22) 


where  the  second  comparison  follows  from  Lemma  5.  Combining  (21)  and  (22) 
readily  leads  to  (19).  ■ 


5  Proposition  4  in  action 


Proposition  4  will  now  be  used  to  discuss  the  bounds  of  Rasmussen  et  al.  [7]  and 
of  Mao  and  Habibi  [3,  4,  5].  Given  N  independent  on-off  sources  (Pn,  Bn(fn )) 
(n  —  1, . . . ,  N ),  all  these  results  express  bounds  of  the  form 

N  L 

PnBn(fn )  <cx  Pnew  ^M/new)  (23) 

n= 1  £=1 


with  i.i.d.  on-off  sources  (Pnew,  -B^(/new))  (£  =  1, . . . .  L)  for  appropriate  con¬ 
stants  Pnew  >  Pm  ax  and  /neW  in  (0,1],  and  positive  integer  L  <  N. 

This  upper  bound  (23)  will  flow  from  Proposition  4  if  P*  and  f*  are  selected 
according  to  (17)-(18),  so  as  to  ensure  the  identification 

NP* 

Pnew  —  — ~jP —  (24) 


and 


fn 


r 


once  L  <  N  is  chosen. 


^total 

NP * 


^total 


TPn 


(25) 


5.1  The  bound  by  Rasmussen  et  al.  [7] 


Rasmussen  et  al.  assume  Px  =  . .  .  =  PN  =:  Pc,  and  arbitrary  activity  factors 
fi, . . . ,  fN,  so  that  ?7i tc,tai  =  G  En=i  fn ■  Applying  Proposition  4  with  P*  =  Pc  = 
Pmax  and  L  =  N,  we  find  that 


'^total 

NP* 


1 

TV 


E/n 

n=l 
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and  the  bound  of  Rasmussen  et  al.  is  obtained  in  the  form  (23)  with  L  =  N, 
/new  =  /*  and  Pnew  =  Pc.  As  should  be  clear  from  Lemma  3,  this  bound  is 
simply  a  well-known  stochastic  comparison  result  for  sums  of  Bernoulli  rvs  due 
to  Hoeffding  [6,  p.  359]. 


5.2  The  first  bound  by  Mao  and  Habibi  [3,  Thm.  3] 

We  are  in  the  homogeneous  case  with  P\  —  ...  —  PN  Pc  and  m  1  =  ...  = 
mN  =:  mc.  Consequently,  we  have  Pmax  =  Pc,  mtotai  =  Nmc  and 


fn 


mn 

Pn 


n  —  1, . . . ,  N. 


Whenever  P*  and  L  are  selected  so  that  P*  >  Pc  and  L  —  for  some  positive 
integer  U,  it  is  plain  that  L  <  N  while  (18)  yields 


'^total 

NP* 


mc 

P*‘ 


Now,  select  P*  >  Pc  so  that  ALI  —  UPC;  this  is  always  feasible  by  taking  P*  = 


U  rN 


| "77I  Pc-  Applying  Proposition  4  under  these  conditions,  we  get  the  upper  bound 


N  I  U 


of  Theorem  3  in  [3]  in  the  form  (23)  with  L  =  [ A],  )new  —  f*  and  Pnew  =  U Pc. 


5.3  The  second  bound  by  Mao  and  Habibi  [4,  Thm.  3] 

In  [3,  Conjecture  1],  Mao  and  Habibi  propose  an  upper  bound  that  combines  their 
earlier  bound  with  that  of  Rasmussen  et  al..  This  second  bound  is  discussed  as 
Theorem  3  in  both  [4,  p.  127]  and  [5],  and  deals  with  N  arbitrary  independent 
on-off  sources  (Pn,  Bn(fn ))  (71  =  1, ,  N). 

Pick  an  arbitrary  target  value  Pnew  >  Pmax-  The  conjectured  bound  is  of  the 
form  (23)  with  L  =  Lmh  where 

Lmh  =  r^l  with  Ptotai  =  Pi  +  . . .  +  Pat  (26) 

t  new 

and 

_p  _  ^total 

J new  T  p  •  \A  • ) 

^MH  Tnew 

At  first,  we  try  to  check  whether  this  second  upper  bound  does  indeed  flow 
from  Proposition  4  by  appropriately  selecting  P*  and  /*  according  to  (17)-(18) 
when  L  =  CMH:  Obviously  LMH  <  N  since  Ptotai  <  NPmax  and  (27)  is  just  (25). 
Next,  the  value  of  P*  which  meets  the  requirement  (24)  obviously  has  to  be  given 
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by  P*  =  ^jr-Pnew  Unfortunaly,  the  constraint  P*  >  Pmax  is  usually  not  satisfied 
as  can  be  seen  upon  rewriting  it  in  equivalent  form  as 


NPmax  . 


^  r  total  “I 

-  I  ~p  r 

-*■  npw 


Pn, 


(28) 


In  fact,  if  (28)  is  to  hold  at  all,  it  can  only  hold  as  an  equality,  in  which  case  we 
must  have  P*  =  Pmax.  Thus,  in  general,  the  conjectured  bound  of  Mao  and  Habibi 
is  not  a  consequence  of  Proposition  4  unless 


PPrnax  r  Ptotal  -> 

_ P  _  I  ~P  I' 

1  new  1  new 


(29) 


5.4  A  new  provably  correct  upper  bound 


The  state  of  affairs  just  uncovered  leads  us  to  suspect  that  the  conjectured  upper 
bound  might  be  in  error,  and  indeed  the  end  of  Part  2b  of  the  proof  given  in  [4,  p. 
137]  appears  to  be  in  error.  As  will  become  apparent  shortly,  the  bound  proposed 
earlier  is  “too  tight”  with  too  few  terms. 

To  see  this,  we  recall  that  the  comparison  (19)  will  lead  to  (23)  with  a  pre¬ 
scribed  target  value  Pnew  >  Pmax  upon  selecting  P*  and  /*  according  to  (17)- 
(18),  and  a  positive  integer  L  <  N,  so  that  (24)  and  (25)  hold.  Thus,  once  Pnew  is 
prescribed,  the  integer  L  <  N  and  the  auxiliary  variable  P*  determine  each  other 
via  the  relation 


L 


NP * 


(30) 


under  the  constraints  Pmax  <  P*  <  Pnew- 

Reductions  in  computations  are  achieved  by  selecting  the  smallest  admissible 
value  of  L,  say  Lmin,  so  that  P*  given  via  (24)  satisfies  (17).  These  constraints 
yield 

-^min  •  min{L  =  1, . . . ,  N  :  —  Pnew  >  Pmax}  (31) 


whence  Lmin  and  the  corresponding  P*  are  now  given  by 


P  T 

T  r  at  1  max  i  1  7-,*  -^min 

Lmin  =  |  N— —  |  and  P  =  — — - 

-L  npw  4  V 


-  P 

1  ni 


(32) 


so  that 


r 


Ttl  total 

NP* 


U^total 


T  P 

in  i  n 1  new 


fn 


(33) 


The  new  upper  bound  is  therefore  of  the  form  (23)  with  Lmin  terms.  We  note  that 

LmU  A  Pmin- 
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